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Abstract 

Some spatially homogeneous Bianchi type I cosmological models filled with 
homogeneous "electric" p-form fields are shown to mimic the never-ending oscil- 
latory behaviour of generic string cosmologies established recently. The validity 
of the "Kasner-free-flights plus collisions-on-potential-walls" picture is also il- 
lustrated in the case of known, non-chaotic, superstring solutions. 



1 Introduction 



It has been recently pointed out that the generic inhomogeneous solution, near a 
cosmological singularity, of the low-energy bosonic field equations of all superstring 
models (D = 10, IIA, IIB, I, hets, hetso), as well as those of M-theory (D = 11 
super gravity), exhibits a never ending oscillatory behaviour of the Belinskii- 
Khalatnikov-Lifshitz (BKL) type @. As a rule, the existence of such an infinitely 
continued oscillatory behaviour delicately depends both on the full field content of 
the considered theory, and on the allowance for generic spatial inhomogeneities. Non 
generic models, obtained, e.g., by truncating the field content, and/or by imposing 
some homogeneity ansatz, may turn out to exhibit a finite number of oscillations, 
ending with a monotonic Kasner-like power-law approach to the cosmological sin- 
gularity. Such solutions are of measure zero within the full set of solutions of the 
complete theory and do not contradict the finding of Ref. [Q] which concerns generic 
solutions. 

The study of specialized cosmological models, in particular spatially homogeneous 
models, has played an important role in theoretical cosmology. Of particular interest 
has been the study of the homogeneous Bianchi type-IX model , [|J , for Einstein's 
vacuum equations in D = 4. This model captures, by means of a simple set of ordinary 
differential equations (ODE), the essential features of the generic, inhomogeneous 
BKL oscillatory behaviour ||. One may then wonder whether there exist also in the 
superstring context, simple, homogeneous models which capture by means of ODEs 
the essential new features of the oscillatory behaviour found in Jl| . The purpose of the 
present paper is to answer this question in the affirmative. The simple homogeneous 
models in question are just the Bianchi type I models with appropriate p-form sources. 

Let us first recall the basic results (and the notation) of Jl|. We consider models 
in spacetime dimension D of the general form 



S = J d D x^g 



We have written (Q) in the Einstein conformal frame, and used a specific normalization 
of the kinetic term of the "dilaton" (p. The integer p > labels the various p-forms 
A p = A^,,,^ present in the theory, with field strength = <9 W ± p 

permutations. The real parameter A p measures the strength of the coupling of the 
dilaton to the p-form A p in the Einstein framef]. In the case of M-theory, the dilaton 
is absent, and one must set (p = in ([!]). 

Following ||, we consider a "generalized Kasner" metric g^ u dx^ dx u = —N 2 (dx°) 2 + 
Z)f=i a i { ujt ) 2 i where d = D — 1 denotes the spatial dimension and where u> l (x) = 
e Z j(x) dxi is a <i-bein, whose time-dependence is neglected compared to that of the lo- 
cal scale factor aj. When working in the proper-time gauge, N dx° = dt, the leading 
Kasner-like approximation to the solution of the field equations for g^ and ip derived 
from ([!]) is, as usual @ In a« ~ Pi(x) In t + b^x), <p ~ p^x) In t + ip(x) . The 

1 Actually, the superstring actions (in D = 10) are slightly more complicated than Eq. (|l]), in 
that they include additional couplings between the form fields (e.g. Yang-Mills couplings for p = 1 
multiplets, Chern-Simons terms, {(IC2 — C0 d i?2) 2 -type terms in type IIB). However, these additional 
terms do not qualitatively modify the generic BKL behaviour discussed in [jlj. Similarly, there is a 
Chern-Simons term in D = 11 SUGRA. 
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spatially dependent Kasner exponents pi (x), p^ (x) must satisfy the famous Kasner 
constraints (modified by the presence of the dilaton): 

pj + £p?-(x>) =0, Eft = i. (2) 

i=l \i=l / i=l 

The set of parameters satisfying Eqs. (H) is topologically a (d— l)-dimensional sphere: 
the "Kasner sphere". When the dilaton is absent, one must simply set p v to zero in 
(0). In that case the dimension of the Kasner sphere is d — 2 = D — 3. 

The approximate generalized Kasner solution is obtained by neglecting in the field 
equations for g^ u and <p: (i) the effect of the spatial derivatives of g^ u and <p, and 
(ii) the contributions of the various p-form fields A p . As in the usual BKL approach 
[see for a summary of the evidence supporting the BKL picture] the approach of 
[[j], U assumes that the Kasner-like solution is approximately valid during a sequence 
of "free flight" evolutions of and <p, which are interrupted by "collisions" against 
some "potential walls" associated to the momentarily growing effect of either (i) the 
spatial derivatives of g^, or (ii) the "electric" or "magnetic" contributions of one of 
the p-form fields A p . 

It was found in |], H that the "potential walls" (conjecturally) responsible for 
the collisions are, in the gauge iV = ^fg, of the form V(j3^) ~ J2 Cj exp(2 wj{j3^))^ 
where Wj{0) = Wj fl j3 tl are some linear forms in f3 l = In a«, (3° = (p. Each elementary 
potential wall Vj(/3 At ) ~ exp(2 Wj(f3 ti )) in the (N = y/g) Hamiltonian constraint 
corresponds to one of the potentially "dangerous terms", when t — > 0, among the 
neglected contributions of the types (i) or (ii) mentioned above, i.e. among t 2 i?*, or 

among t 2 T® A ^ and t 2 T[ A ^. Here, ^ denotes the d-dimensional Ricci tensor whereas 
T^ v denotes the stress-energy tensor of the p-form. More precisely, upon replacing 

the "ingoing" approximate Kasner behaviour, we have Vj(/3 M ) ~ t 2wj ( p ) where the 
exponent of t 2 is simply wj(p) = wjop v + wjiPi. This result exhibits the simple link 
between the dominant elementary potential walls Vj{(3^) ~ exp(2wj /J /3 M ) responsible 
for the collisions, and the set of "stability exponents" wj{p) = wjop^+wjiPi discussed 
in 0. 

It has been known for a while ||, ]7) that spatial inhomogeneities in g^ v give rise 
to the following set of gravitational exponents g^ip) = 2pi + J2e^i,j.k Pe / j, «' 7^ k, 
j 7^ k). Ref. [|I|] found that the presence of p- forms gave rise (for each type of p-form) 
to two additional sets of potential walls and stability exponents: the electric exponents 
e< h...i p (p) = Ph +Pi 2 + ' ' ' + Pi P ~ \ \Pip (where all the indices i n are different) and the 

magnetic exponents j d 1 (p) = pj l + Pj 2 + ■ ■ ■ + Pj d ^ p _ 1 + \ \ p^ (with all indices 

j n being different). In the case of the dilaton- free models, one must set p^ to zero in 
these expressions. 

The essential new result of |l| was the finding that, in all low-energy super- 
string models (D = 10, IIA, IIB, het E , het so ; and D = 11 SUGRA) there ex- 
ists no open region of the Kasner sphere (fj) where all the exponents {wj(p)} = 
{9ijk(p),^...i p (p)^ ( f^j d _ p _ 1 (p)} can be simultaneously strictly positive. [We shall 
give below a simple direct proof of this fact for D — 11 SUGRA.] In other words, 
the Kasner-like evolution cannot stay monotonic, but will always be "deflected" by 
at least one potential wall to a new Kasner regime characterized by new Kasner ex- 
ponents p'i(x), p' v (x) (given by the universal collision law of [|TJ) as well as, in general, 
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new "Kasner axes" e'Hx). This new Kasner solution will again be deflected because 
there is always one "dangerous" potential term that takes over, since at least one wj(p) 
is negative, etc. This never ending oscillatory behaviour is similar (though more com- 
plex because of the greater variety of potential walls) to the BKL oscillations in D = 4 
pure gravity, but contrasts very much with the ultimate monotonic Kasner behaviour 
holding for pure gravity in D > 11 0, and for the Einstein-dilaton system in any D 
||. As pointed out in |]J it is the presence of various form fields (e.g. the three form 
in D — 11 SUGRA) which provides the crucial source of generic oscillatory behaviour. 
This leads us to propose in this paper to consider simple homogeneous models which 
contain no "gravitational walls" but which contain sufficiently many "form walls" to 
prevent the possibility of an ultimate monotonic Kasner behaviour. We view these 
models as analogs of the Bianchi IX model, i.e. as toy ODE models for studying in 
detail the chaos induced by the form fields. 

2 M-theory and electric potential walls 

We first consider the case of M-theory. In that instance, the toy ODE models 
which exhibit the required features are in fact the simplest "Bianchi I" type ver- 
sions of the action (|]). These models contain: (i) no dilaton, (ii) a generic spa- 
tially flat .D-dimensional metric with closed d-bein forms duo 1 = 0, or equivalently 
ds 2 = —N 2 (x°) (dx ) 2 + gij(x°) dx l dx^ = 1, . . . , d = D — 1); and (iii) a generic, 
homogeneous p-form potential A = (1/p!) A^.; (x°) dx %1 A ... A dx lp . [In M-theory, 
p — 3, but it is of interest to leave p unspecified at this stage.] 

These models have several useful features: (a) because the potential is homoge- 
neous, the curvature (p+ l)-form F = dA = do A^...^ dx° A dx 11 A ... A dx %v is purely 
electric; there is no magnetic field; (b) the vanishing of the structure constants C l - k 
of the homogeneity group allows the purely electric (p + l)-form F to automatically 
satisfy the Gauss constraint (when Cj fc ^ 0, see |J); and (c) when D = 11 and p = 3 
this model is equivalent to (Bianchi I) SUGRA because the vanishing of the magnetic 
part of F means that F A F vanishes, so that the addition of the Chern-Simons term 
A A F A F to (|l|) has no effect in the field equations. 

After discarding a total derivative and a trivial volume factor, the action ([!]) 
(reduced by the ansatze (i),(ii) and (iii)) reads, in Hamiltonian form 

S H = / dx°(^ gij + -^K...^ ~ NH) , (3) 
J p. 

where ir l J and 7r^"' ip are respectively the momenta canonically conjugate to gr^ and 
Ai 1 ...j , and where is the rescaled lapse of weight minus one related to the standard 
scalar lapse through N = N/ yfg. In (^), TC is explicitly given by (^ = gj k 7i z g k ) 

n = 4, - «) 2 • ^ g hjl . . . g ipjp 7#"* . (4) 

The dynamics of Bianchi I models admits many first integrals. Indeed the action 
S is invariant under the following two sets of rigid symmetries: (i) an arbitrary 
GL(d) transformation, described by a (non symmetric) d x d matrix Aj acting by: 
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9ij = K^ b j9ab, A' h ip = Ajl . . .Aj A h „. jp , and (ii) an arbitrary shift of the p-form: 
^h...i p = ^h-ip + The corresponding Ncether conserved quantities are 

-pi _ o-i I _ isi...s p -i A ch...ip _ h-ip 

H 3 - Z7r 93 + (p-l)\ A t A \P) 

[In Eq. d^) and the following, we reserve the notation £••• to denote a numerically 
fixed constant of motion, while ir^ denotes a conjugate momentum.] By inserting 
these integration constants into the definition of the momenta in terms of the first- 
order time derivatives of g^ and A^...^, one gets a first-order system in (g, A) of the 
symbolic form d T g = P\^(g, A), d T A = P p (g), where Pi^(x, y) is a polynomial of first 
degree in x, and first degree in y, and P p (x) a polynomial of degree p in x. 

In addition to the conserved quantities (|), the dynamics is subject to the con- 
straint 7i = which one obtains by varying N in the variational principle. This 
constraint is preserved in time and defines a further constant of the motion, which is 
generically independent from the previous ones. 

Rather than investigating the polynomial equations of motion obtained by elim- 
inating the momenta, it is more instructive for our purposes to analyse the reduced 
Hamiltonian for the dynamics of the metric variables, after elimination of the A's 
(which are "ignorable" coordinates) and given the constants of motion £. Actually, 
it suffices to replace by £ •" in the Hamiltonian (B). 

We work as usual in the r-gauge, defined by N = 1, i.e., dr = dt/^/g. We 
introduce the positive-definite "potential" for (given the constant "electric field 
strengths" £ h - A p) 

v £ ^\g) = l^gi^---9ip jp £ h -- Ap £ jl --- jp - (6) 

z p. 

Then the dynamics for g^ follows from the r-time Hamiltonian 

H?> (g, k 9 ) = (g ia g jb - 9a g ab ) < *f + v£ p) ( 9ij ) . (7) 

The initial conditions must be chosen such that the zero-energy condition H^' = 
is satisfied initially (it is then preserved by the evolution). Actually, the reduction 
leading to (0) is far from optimal because we have not taken into account all the 
symmetries of the problem, i.e. all the constants of motion, but it is convenient for 
our present purpose which is to give simple examples of the analysis of [Q. 

Let us qualitatively discuss, a la BKL ||, or rather (as we work in a Hamiltonian 
framework) a la Misner ||, the dynamics generated by (|7]). I.e., let us verify the con- 
sistency of a picture in which gij{r) evolves by a succession of "free flights" interrupted 
by "collisions" on the potential wall Vr P \g). The free flights are the periods where 
V^\g) is much smaller than the kinetic energy terms ~ g 2 ~ 7r| in (|7|). When this 
is the case, gijij) undergoes a simple geodesic flow in the supermetric defined by the 
Hamiltonian Hf \ i.e. (0) without the potential term. The invariance of under 
GL(d) gives, as above, the constants of motion VL = 2 7t^. The free flight evolution 
of g is therefore given by d T gij = g ia Vgj — (d— 1) gij Vg a . This is a linear equation in 
gij. It is easily solved by diagonalizing the constant matrix V^, say Vgj = V l g <5* after 
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some GL(g?) transformation V g — > A x V g t±. The solution for gf^ (in the new linear 
frame) reads c/g eeflight (r) = exp (2/3j 0) + 2v < r) with 2u* = 7>j - (d- l)" 1 ^P*. 

In terms of the cosmological time t = J y/gdr, this is a Kasner solution with 
a,i oc t Pi and Kasner exponents equal to p, = v l / i^2, s v s ). This (approximate) Kasner 
solution can continue, uninterrupted, if and only if all the terms in the potential 
V<P*(g) tend to zero as t — > 0. [Note that a typical contribution to the kinetic energy 

terms ~ (d T a,i/ai) 2 is of order unity, i.e. independent of r as r — > oo, i.e. t — > 0.] 

j 

In the special diagonal basis where we have written the solution, the electric field 
gn...ip w jjj h ave j n generic solutions (i.e. excluding special solutions of zero measure), 
non zero entries for all its components. Therefore the potential ([]) is a sum of positive 
terms of the form a? . . . a 2 (£ J1 - 4 p) 2 ^ ^(p«i +-+«,,)_ The Kasner stability condition, 
i.e. the fact that all such terms tend to zero with t, is therefore precisely the condition 
pointed out in |] that the electric exponents (with p v = in the present dilaton 
free model) be strictly positive for all possible choices of (different) indices ii, . . . , i p . 
This condition is equivalent to e^ in (p) > for some p on the Kasner sphere, where 

e min ip) = Pi + P2 + • • • + Pp j with pi < p 2 < . . . < Pd, denotes the smallest electric 
exponent associated with the presence of a p-form. 

3 Electric chaos and spacetime dimension 

In the pure electric Bianchi I toy model, we have lost all the other stability conditions 
linked to spatial inhomogeneities in the gravitational field, or to the generic presence 
of magnetic-type field strengths. Let us, however, delineate for which values of the 
spacetime dimension D = d + 1, and of the degree p of the form the electric stability 
conditions (without the dilaton term) are sufficient, by themselves, to imply a chaotic 
behaviour, in the sense of a never ending oscillatory behaviour, as t — > 0. We shall 
prove the following (we consider only d > 3, i.e. D > 4, so that the Kasner sphere 
exist as a continuous manifold, and forms of degree p < 3): 

Theorem 1. Let p\ < P2 < • • • < Pd be ordered Kasner exponents running over 
the Kasner sphere S d ~ 2 (Y,Pi = Hp 2 — 1)> (with d > 3), and let e^ in (p) denote the 
smallest electric exponent associated to the presence of a p-form: 

(i) in the case of a one-form (p = 1): = pi can never be > on S d ~ 2 ("electric 
chaos" for any D > 4), 

(ii) in the case of a two-form (p = 2): if D = 4, ej^ = P\ + P2 can be > on S d ~ 2 
but if D > 5, e m ( n can never be > ("electric chaos" in D > 5), 

(iii) in the case of a three-form (p = 3): if D < 6, e ffi n = pi + p 2 + p 3 can be > on 
S d ~ 2 , but if D > 7, e^; n can never be > ("electric chaos" in D > 7). 

Note that a consequence of (iii) for D = 1 1 is the result announced in [jlj , namely 
the fact that the 3-form of SUGRA creates, by its sole electric effect, chaos in D = 11 
supergravity. [We sketched in QXJ another proof (which was our original proof) of that 
fact. The proof we give below is much simpler.] Note also that the chaotic nature of 
a one- form in D = 4 (and 5) was clearly recognized in Ref. [|i~0[|. 
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We shall only prove the more difficult part of the theorem, namely, the assertion 
(hi). The other claims are proved similarly ||. First, we note that: (a) if d = 3, 
e min = Pi + P2 + P3 = 1 is always > on S 1 ; (b) if d — 4, e$ m = p% + p 2 + 
P3 = 1 — is > almost everywhere on S 2 ; and (c) if d — 5, the particular point 
(p( 0) ) = (-l^l^lj on S* 3 satisfies e^ n (p^) = | > 0, so that there exists an 

open neighbourhood of p^ where e^ m (p) > 0. Let us now prove that, when d > 6, 
the assumption e^; n = pi + p 2 + P3 > leads to a contradiction. With our convention 
of ordered p,'s, this assumption implies < p 3 < p 4 < . . . < pd- It follows from the 
Kasner conditions that the double sum K{p) = J2i<i<j<d PiPj — vanishes. Let us 
distinguish in this double sum the indices i, j = 1 or 2 from the indices > 3, which we 
denote by a, (3 = 3, 4, . . . , d. [The indices a, (3 take d — 2 values.] After some simple 
rearrangements, one easily checks that the following algebraic identity holds: 



K{p) = G» 3 - Pl)(P3 - Pi) + (Pi +P2+ P3) 
P3+P3+ Y] Pa 



P3+ 

3<a<d 



P3 



3<a<d 



3<a<0<d 



The first term on the R.H.S. of (||) is clearly > because of the ordering of the p^s. 

(3) 

The second term is strictly positive if the assumption > holds somewhere. 
Now, if the number of terms ((d — 2)(d — 3)/2) in the last double sum is larger or 
equal to the number of terms (2 + d — 2 — d) in the (expanded) penultimate simple 
sum, the inequalities < p% < p^ . . . < pd show that the difference between the two 
will be > 0. This would lead to the contradiction K(p) > 0. Therefore > is 
contradictory when (d — 2){d — 3)/2 > d, i.e. (as easily checked) d > 6. Q.E.D. 

Theorem 1 tells us which simple "electric Bianchi I" models incorporate enough 
potential walls to mimic, within a simplified context, the form-induced chaos that 
Ref. found to hold in all (inhomogeneous) superstring models. The chaotic ho- 
mogeneous model relevant to M-theory is the dilaton-free Bianchi I model with a 
3- form (p = 3), in spacetime dimension D = d + 1 = 11. The study of this model, 
with a generic electric field, might teach us something about the nature of the chaos 
induced by the 3- form in SUGRA n . We emphasize that it is essential to consider 
a generic solution to have chaos. There exist particular solutions (of measure zero 
among all solutions) which exhibit only a finite number of oscillations. For instance, 
solutions with just one collision are given in 0. Theorem 1 shows that a purely elec- 
tric, homogeneous 3-form leads to chaotic oscillations in any space dimension d > 6. 
Therefore a study of this simple model in d — 6 (which has many less variables than 
its d = 10 analog) might be sufficient to learn about the nature of the chaos induced 
by the 3-form. Actually, Theorem 1 allows us to consider even simpler models (with 
less dynamical variables) by lowering both the degree of the form and the dimen- 
sion. The simplest models that one can consider to study the nature of the chaos 
induced by a form field (and to compare, and/or contrast it, with the Bianchi IX 
chaos, which is induced by the gravitational field) is a vector field p = 1, i.e. a usual, 



Maxwellian electric field, in a (^-dimensional Bianchi I model with d > 3 [10|. The 



case d = 3 has been studied in |Il], |L|, [I3[ with the conclusion that, indeed, a generic 



homogeneous electromagnetic field induces chaos. Note, however, that in d — 3 the 
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electric and gravitational exponents characterizing the walls are proportional so that 
they lead to identical collision laws. It would therefore be more interesting (to un- 
derstand the specificities of the form-induced chaos, as well as the effect of increasing 
the dimension) to study the case d > 4. 

4 String theories and chaotic cosmological models 

From the above findings for D=ll SUGRA, one can easily construct homogeneous 
cosmological models with an infinite number of Kasner oscillations for type IIA string 
theory. The above models can, indeed, be interpreted as Bianchi I models in ten 
dimensions, with a general metric, a Kaluza-Klein vector field, a dilaton, as well as 
homogeneous 2-form and 3-form fields. Thus, Bianchi I models with homogeneous 
potentials (i.e., only electric fields) exhibit an infinite number of oscillations in type 
IIA string theory. By T-duality, the results extend to type IIB, but since some of 
the magnetic fields now do not vanish, there is no action principle with homogeneous 
p-form potentials from which the equations of motion derive (some of the potentials 
must be inhomogeneous) . 

One can also devise Bianchi type I models which exhibit a never-ending set of 
oscillations for the (D = 10) heterotic and type I string theories because the gravi- 
tational walls related to spatial inhomogeneities (which are absent in the Bianchi I 
context) are not necessary to induce the oscillations. As shown in 0, the electric 
stability conditions associated with the 1-forms and the magnetic stability conditions 
associated with the 2-form cannot be simultaneously fulfilled in D = 10. So, by con- 
sidering homogeneous sources of this type, one gets cosmological models that exhibit 
the required behaviour^. 



5 A non-chaotic model 

The dilaton and the p-forms (p > 0) have quite different effects on the oscillatory 
behaviour of cosmological models. While the dilaton modifies the Kasner exponents 
in a way that they can all be positive - and so, tends to stabilize the Kasner behaviour 
- the p-forms (p > 0) do not change the Kasner relations but induce collisions - and 
so, tend to destabilize it. This observation is useful even for solutions with a finite 
number of oscillations (which form a zero-measure set) as it enables one to understand 
the dynamics of these models in terms of free-flight motions interrupted by collisions. 
We illustrate this feature with a discussion, for comparison purposes, of the case 
of a two- form (p = 2), in presence of a dilaton (coupled in the bosonic-string or, 
equivalently, heterotic-string way). This case is interesting because it is known to 
be exactly integrable We wish to show, on this example, how its dynamics can 



be described in terms of BKL-like collisions, and how this collision analysis correctly 
predicts: (i) the presence of only a finite number of oscillations, and (ii) what is the 
final "out state" after the oscillations have ceased. The model we consider here is 
described by the action ([!]), in any spacetime dimension D = d + 1, with p = 2 and 

2 Magnetically induced chaos was pointed out in a particular cosmological context in JTJJ]. It was 
also pointed out there that the dilaton makes this magnetic chaos disappear. However, if one brings 
in also the electric fields of the 1-forms, one gets chaos back again. 
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Afp=2) = —4/y/D — 2. This specific value of the dilaton coupling in the Einstein frame 
follows from starting from a string frame action of the form 



S = J d D x ^ e~* \r(G) + (V0) 2 - ^ H Xllv H x ^ 



(9) 



where = exp(20/(D — 2))g^ u is the string frame metric, = y/D — 2 the 
dilaton (normalized such that the string coupling squared g 2 s = e^) , and where = 
d\ B^ u + B u \ + d u B\fj, is the exterior derivative of the 2-form B^ u . The action (|9]) 
is (in the corresponding critical dimension) a common subsector (NS-NS) of all string 
theories. [But it is only for the closed bosonic string (D c = 26) that (§) is the full, 
tree-level bosonic action.] Anyway, we consider here @ as a toy model to see how 
dilaton couplings can modify the stability exponents of the p-forms so much as to 
quench the form-induced chaos. To study the Kasner stability of the model @ it 
is convenient to work with the string- frame Kasner exponents, say ai, i = 1, . . . , d, 
instead of the Einstein-frame ones p^. The string-frame exponents are defined by 
writing the Kasner-solution in terms of the string-frame scale factors aj and the 
string-frame cosmological time t {G^ v dx^ dx u = —dt 2 + J2 (fli dx 1 ) 2 ^ : a, oc t a \ They 

are linked to the p^s by pi = ((d— 1) a% — a)j{d — 1 — a), = [yd — la)/(d — 1 — a), 
where a = (Ya a i) ~ 1- m terms of the a's the Kasner sphere S d ~ l defined by Eqs. (|2|) 
becomes simply Y^=\ a i — 1- It is easily found, either by transforming the Einstein- 
frame exponents, or by a direct analysis in the string frame, that the string-frame 
exponents e, defined such that the "dangerous terms", i.e. the potential walls, grow 
oc t 2e , (e = (1 — a/(d— l))e, for each corresponding Einstein-frame exponent, e) 

12) — (2) 

read g^ k = 1 + ai — otj — e\j = «j + acj and j = 1 — ai — otj — a^- Here, 
ijk£\ . . . id-3 is a permutation of 1 ... d. We shall show in that, when considering 

all these exponents, the Kasner stability conditions (~g,~&( 2 \ b all > 0) can be satisfied 
only when D > 11. This implies, for instance, that the NS-NS sector of type II or 
heterotic superstring theories in D = 10 are chaotic. However, in the present paper 
we are interested in considering only the simple purely electric Bianchi I models. In 
these models, there are only "electric walls", so that the only stability condition to 
satisfy is e-^ = ai + ctj > 0. Clearly this is always satisfied in some region of the 
Kasner sphere X)f=i a f = 1 ( as long as d > 2, which is anyway necessary to consider 
a Bij). Therefore, we conclude that the electric Bianchi I model ([5]) (in d > 2) will 
not be chaotic, and will contain at most a finite number of oscillations. 

We can, in fact, be more precise and determine the maximum number of os- 
cillations by considering the "collision law" induced by an electric wall. This law 
follows from the general collision law given in |]J (Eq. (9) there). The string- frame 
Kasner exponents are the "velocities" a« = djT/df, where j3 = lna^, dr = Xe^di, 
with a normalization constant A such that ao = d<j)/dr = — 1. In terms of the 
incoming (aj) and outgoing (a-) exponents, the collision law corresponding to the 
wall oc e 2 ^ 1+/32 ^ oc e 2 ^ ai+(X2 ' T (see below) reads a[ = —a 2 , a' 2 = — a±, a' a = a a 
(a = 3, 4, . . . , d). If one starts at some initial time r with some initial values of the 
Sj and with some electric components £ ^ that are all of the same order of magnitude, 
one generically expects that the first collision encountered as r decreases^ will be the 

3 We follow the evolution toward the cosmological singularity at t = 0, say starting from t < 0, 
with t r~j —<f> ~ + In |f | going to — oo at the singularity. 



s 



one associated to the fastest growing wall (if there are such walls), i.e. to the most 
negative a, + atj 7 because afa 2 - oc e 2 ^ i+ ^ oc e 2 ( ct *+ a j) r . We recover the notion of 

(2) 

electric 2-form stability exponent e^- = oti + ctj. Let us order the Kasner exponents 
as «i < «2 < • . . < aa- We expect the following generic qualitative evolution: (i) if 
the smallest e\j , i.e. a,\ + 0^2, is < a collision will occur against the growing wall 
oc a^af, and the outcome of this collision will be the unordered set a\ given by the 
above collision law. Note that a[ + a' 2 = — (a% + 0:2) has become > 0. However, they 
may remain some < electric exponents associated to the unchanged a' a = a a \ then 

(2Y 

(ii) if the next smallest e\j , namely 0:3 + 0:4 is < 0, a second collision will occur with 
the effect of reversing the signs of both 0:3 and 0:4. This process will continue until all 
= cti + oij have become > 0, which means that at most a° nt can be < 0, all the 
other ones (here a° ut < o;^* < . . .) ending up being > 0. If we do the same reasoning, 
starting from Tq, but going in the sense of increasing r (i.e. towards \t\ = 00) we 
conclude (mutatis mutandis) that, near \t\ = 00, the "incoming" values a- n must all 
be < 0, except maybe a™ which might be > 0. As each collision changes the sign of 
two Kasner exponents, the total number of collisions cannot exceed [d/2], where [. . .] 
denotes the integer part. 

The qualitative picture just explained can be confirmed by an exact analysis. 
Indeed, it happens that the model at hand is exactly integrable []15| |. However, the 



form of the "general solution" given in Refs. |15|| is not optimally useful because it 



contains a constant (2d) x (2d) matrix A which must satisfy a whole set of non linear 
constraints. We found, however, a simpler way of writing the general solution. It is 
enough to use two facts: (i) a particular, exact solution containing 2d — 1 arbitrary 
constants is known, namely the Kasner solution 

G%(r) = exp(2^ (0) + 2a t r) 5 tJ , B% = ; (10) 

and, (ii) there is an 0(d,d) symmetry transforming solutions into other solutions. 
Using a (d x d) block notation for (2d) x (2d) matrices the 0(d, d) group is realized as 



U — ( y g j, constrained by U 77 U = 77, where 77 = ( ^ ^ J . This is equivalent 

to the constraints W T Y + Y T W = 0, X T Z + Z T X = and W T Z + Y T X = 1. 
The action of U G 0(d,d) on the matrices G = (G. t j) and B = (Bij) is defined 
by M ncw (r) = U M old (r) U T , where the symmetric (2d) x (2d) matrices M are con- 
structed as ] 

M=(!?~ -'TJLV fir 



By looking at the action of infinitesimal transformations U — 1 + u, with u 

w x 



KBG- 1 G-BG- l B 

finitesimal transformat: 

rp J , with x and y antisymmetric, and w arbitrary, one finds that: w induces 

a trivial global change of linear frame, y induces a rather trivial shift of B^ by con- 
stants, while x generates (when acting on ([II])) a time-dependent 5 Bij(r) of enough 
generality to match any initial data 5Bij(T Q ). Going back to finite transformations 

/ 1 X 

U, we conclude that it suffices to act by a triangular 0(d,d) matrix U 



.0 : 

with a generic, antisymmetric matrix X, on the generic time- dependent particular 
solution (|10"D , to generate a solution of the system with sufficient, physical generality. 
Thus, we conclude that the general solution can be written in terms of the Kasner 
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one P!D| ) as 



G 



-i 



Gq — X Gq X , 



(12) 



B 



G X G 



— — <J"0 



Cm X G 



Gq X Gq ( 1 — -AT (jr o X Gq) 



(13) 



We shall focus on the behaviour of the metric, using the explicit form of fll2"|), i.e. 
(with X ij = -X ji ) 




multi-collision picture explained above. However we have verified that (0) is fully 
compatible with the general collision picture, and does indeed contain up to [d/2] 
collisions, which can, for a general set of initial conditions, be well separated from 
each other. First, it is instructive to write down explicitly fll~4] ) in the case where the 
matrix X couples only to a 2 x 2 subblock: G = diag(aQ,6 ) with a (r) oc e aiT , 

2 

&oCr) oc e a2T . With X 12 = x one finds G = diag(a 2 ,6 ) with the new scale factors 
a 2 = a^/il + x 2 alb ), b = b /(l + x 2 alb ). If a^b^ grows toward the singularity 

2 

(i.e. if it corresponds to a growing potential wall Vg oc (£ 12 ) 2 aQ& ) we see that these 
equations describe a collision which turns (a 111 ) 2 ~ a 2 ,, (b ) 2 ~ 6 , into (a° ut ) 2 ~ 
x~ 2 b Q , (b ) 2 ~ x~ 2 a,Q 2 . This corresponds precisely to the collision law given above 
for the (1, 2) block. One can check that this "elementary" collision will approximately 
take place each time a 2 x 2 subblock of Go separates itself from the other ones by an 

2 

especially strong growth of Oq b . We have also done some numerical experiments (for 
d = 4) with a generic, random antisymmetric matrix X and some initial values of the 
Kasner exponents «j ensuring that two of them are especially negative, and the two 
others less negative, and we have found that, in such a case, one indeed witnesses 
two successive collisions, each one reversing the sign of a pair of Kasner exponents. 
[We define some "instantaneous" Kasner exponents as the eigenvalues of the matrix 
K = \ G _1 &f G.} One should, however, say that when there happens to be no clear 
hierarchy between the values of the a's, i.e. no clearly dominant wall among the 
total potential Vs(G), the evolution of the eigenvalues of K looks rather like a single, 
complex collision whose outcome is to change the sign of all dangerous a's at once. 

In this respect, we note that it is easy to find analytically, from ([14]), the net overall 
(S'-matrix) effect of the entire collision process (independently of whether it can be 
viewed as made of several separate, intermediate collisions). Let us first consider the 
case where d is even and where one chooses all the a's in (0) to be negative, i.e. 
where all the clj = grow toward the singularity. It is easy to see from (|i~4|) that 
the incoming state (\t\ — > oo, f — » +oo) is given by G™ ~ Gq^ = (aoi) 2 5ij, while 
the outgoing state (\t\ — > 0, f — > — oo) is generically (using detX ^ generically 
in even d) equivalent (after rediagonalization) to G° ut ~ (<%)~ 2 This shows that 
a° nt = — aj n , in keeping with the result predicted by the "collision" analysis. The more 
general case where d might be odd and where the ordered seed a's, «!<...< a^, 
might be of both signs can also be seen to confirm our collision analysis. 
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6 Conclusions 



In this letter, we have shown that some spatially homogeneous Bianchi type I models 
coupled to appropriate p-form fields are sufficiently complex to mimic the never-ending 
oscillatory behaviour exhibited by generic string cosmologies [|J. This observation 
should open the door to numerical investigations of the effect of the collisions induced 
by the p-forms, which are different from the gravitational ones. We have also shown 
that the picture of free-flight motions interrupted by collisions is useful even for models 
with a finite number of oscillations, as are the solutions of . 
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